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1. Introduction 

In a recent paper pQ , we presented a new class of twisting type N vacuum solutions of the 
Einstein equations with nonzero cosmological constant A. These type N solutions admit 
a twisting congruence of shearfree and null geodesies aligned with the unique quadruple 
principal null direction. They were shown subject to a rather simple-looking second- 
order nonlinear ODE as imposed by the field equations for type N. Various special and 
series solutions were found or constructed from this ODE. In this paper, we move on to 
discuss their physical meanings and show that these new exact solutions can serve as 
models for the behavior of gravitational waves in cosmology. 

Certain aspects of these solutions, such as the conformal factor, conformal infinities, 
periodicity, etc., can be found in the general discussion on algebraically special twisting 
spacetimes by Hill and Nurowski [2] . Here we focus on the local interpretation associated 
with the equation of geodesic deviation. Similar analysis has been done for non-twisting 
type N solutions with A, i.e., Kundt class and Robinson- Trautman class, by Bicak 
and Podolsky [3] (see also [HE]). However, a complication for our solutions is that 
the nonzero twist is associated with certain non-integrability, and therefore causes a 
lack of 2-dimensional wave surfaces in the spacetimes that always exist in non-twisting 
solutions. Hence due to extra cross terms in the metric, the coordinate basis adopted in 
[3] cannot be directly used in our more complicated twisting solutions, and we decide that 
it would be much more convenient to use non-holonomic bases instead, and therefore we 
generalize all derivations to such bases. In fact, we find that all major results involving 
the equations of geodesic deviation remain the same for this generalization. 

In the next section, we review the solutions discovered in [1] and present their Levi- 
Civita connection 1-forms and the only non-vanishing Weyl scalar, \1/ 4. Then in Section 
3, we write down the general geodesic equations in the non-holonomic basis and use 
Killing symmetries to simplify them. In Sections 4 and 5, a frame for a physical observer 
and an associated null tetrad basis are constructed along arbitrary timelike geodesies, 
with respect to which the relative motion of test particles is studied. Conditions are given 
to determine those geometrically privileged geodesies along which the observer's frame 
can be parallel-transported. Hence in Section 6, these privileged geodesies are identified 
as the radial geodesies and determined explicitly. Their properties are discussed in 
Section 7, which are all pointing to a Killing horizon. In Section 8, we calculate 
wave amplitudes along radial timelike geodesies and demonstrate that the observer sees 
gravitational waves decaying exponentially fast, which agrees with the cosmic no-hair 
conjecture. At the end, we make comments on the appearance of A as a proportionality 
factor in the waves amplitudes, which we find quite unusual. 
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2. The twisting type N vacuum spacetimes with nonzero A 

With the real coordinate system x a = (x, J, u, r), the metric found in [I] can be written 

as 

g = 2 (u>W + u, 3 ^ 4 ) , (1) 
where the null tetrad is given by 

w 1 = Rd(, u 2 = Rd(, 

u 3 = R\, uj 4 = R (dr + Wd( + Wd( + HX) 
with real-valued P = P(J) > and P' = dP/dJ such that 

dC = dx + idJ, A du + 2Mx 



(2) 



(3) 



P ' -PdjL 
L=-e~*' /^exp(/^|^dj)dJ, 

>P 

R = 7^ 77 -7T < r < 7T, 

2 cos 2 

W= - ( -P' + AJ + iCi ] (e" ir + l), H = --APcosr. 
2 \2 / 6 

Here C\ is an arbitrary real parameter and the function P(J) must satisfy 

p* _ ( p/ + 2A ^) 2 _ ?C? _ 10 

2P P 3 ' lJ 

as required by the vacuum Einstein field equations with real arbitrary constant A. By 
the transformation J = w/A, P = g(w)/A > 0, with A / 0, the above equation can be 
put into a A-independent form 

„ (g' + 2wf 2CI 10 

9 = 7, ( 5 ) 

2^ g 3 

These ODEs allow solutions for either A > or A < 0, though the construction of the 
explicit general solution, if possible, still remains an open problem. (See [1] for examples 
of special solutions and series solutions.) 

Now we introduce the non-holonomic null basis {e Q } = {ei, e2, e3, 64} dual to the 
tetrad Q2J): 

ei = ^ (9 - Wd r ), e 2 = 1 (8 - Wd r ), 

f R x (6) 

e 3 = -= (do - Hd r ), e 4 = — d r , 

where we define 

d = d ( -Ld u , d ( = ^(d x -iPdj), 

d = i(BL - dL)d u = -P(djL)d u . 

In particular, the vector field e 4 is tangent to a twisting congruence of shearfree null 
geodesies, and is also aligned with the quadruple principal null direction of the metric. 
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To ensure a nonzero twist along this congruence (also for the basis {e a } to be valid), it 
is required that 

\{BL - dl) = -P(djL) y£ 0. (8) 

With all said, we present components of the Levi-Civita connection 1-forms r A = 
r A (jJ v calculated from Cartan's structure equations &uj x + r A ^ A oj m = for the null 
tetrad (ED: 

riai = -r 122 = ^ [(e" ir + l)P' + 2(e" ir - 1)(AJ + id 

Tm = -^R AP ( 2 + COSr )' = "2^, 

r 231 = —AP ( ie~ ir - tan - + 2i 
231 12R V 2 



r 233 = ^ Ap ( e ~ ir + + 2AJ - 2iC< i)' 



r 



241 



\e- ir ' 2 
2Pcos§ 



(9) 



r 34 i = r 342 = ~ [(e~ ir - 1)P' + 2(e~ [r + 1)(AJ + id)] 

OiX 



343 



r 1 r 

jAP(2 + cosr) tan -, r 3 44 = — — tan- 



12R 2' 2P 2' 

r 2 32 = r 23 4 = r 2 4 2 = r 2 4 3 = r 2 44 = o. 

Those unlisted components can be obtained by either T^ u = — r„ M or complex 
conjugation on (J9j) which interchanges the indices 1 -H- 2 and leaves 3 and 4 unchanged, 
e.g., T 214 = T 124 , T 211 = r 122 . The only non- vanishing Weyl scalar is 

# 4 = C3232 = ~ [AJP' - |AP + 2A 2 J 2 - AC 2 - 2id {P' + 3AJ)] e" ir / 2 cos 3 J 

A ( 10 ) 

= -~ H -%g+ 2w 2 - 4d 2 - 2id (</ + 3w)] e^ 2 cos 3 1 

which, being proportional to A, requires A ^ for type N solutions. In addition, the 
metric has the following two Killing vectors [6] 

d UJ d x -C 1 ud u , (11) 

which we will use to simplify the geodesic equations. 

3. Geodesic equations 

Using the non-holonomic basis ([SI) , we consider a freely falling test particle (observer) 
with the four-velocity 

U a = (u l , M 2 , M 3 , M 4 ), U 1 =U 2 , 

along an arbitrary timelike geodesic such that 

u- u = 2(uV + A 4 ) = e (12) 
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with e = — 1 (also, e = if one considers null geodesies). From the fact that 
u = u a e a = x^dxn as expressed in non-holonomic and coordinate bases, we obtain 

u = -±-[L(u 1 +u 2 )+P(d J L)u 3 }, (13) 

r = — - (Wu 1 + Wu 2 + Hu 3 - u 4 ). 
R 

with ■ = d/dr and r the proper time. The geodesic equations read 
d?/ 1 d?/ 2 

o = — + rvuv, o = — + v x ^u\ 

du 3 „ ,. dii 4 

— + r^u", o = — 



o = — + r^u", o = — + T 3fll/ uV, 



(14) 



from which one can verify, using (151), that ^(m 1 ?/ 2 + m 3 m 4 ) = (hence, the length (fT2l) 
is constant along geodesies). For a Killing vector the product is a conserved 
quantity along geodesies with the four- velocity u^. Hence from (II ip . we find 

ft = z^(^"^). (15) 

C 3 = i? [u 1 + u 2 + + H> 3 ] + 2C 2 L - C 2 C iU {t), 

where C 2 and C3 are real constants. One can show that the above two simpler equations 
are indeed first integrals of the system (114j) . However, with an implicit P(J) generally 
satisfying (T4J), these geodesic equations are still very complicated to solve directly for 
x M (r) without additional assumptions. 

4. Geodesic deviation 

The equation of geodesic deviation reads 

R\ d u^u\ (16) 



where u = dx/dr = u a e a , u u = — 1 as introduced before, and Z(r) is the displacement 
vector. Note that (fT5l) is still valid despite the use of non-holonomic bases [7] . Following 
the construction^] in [3], we set up the observer's frame {e( Q )} along the geodesic with 
e( 4 ) = u and spacelike orthonormal vectors {e^), e( 2 ), e( 3 )} in the local hypersurface 
orthogonal to u, i.e., e (a) ■ e {/3) = fl^e^e^ = r](a)W) = diag(l, 1, 1, -1). The dual basis 
is e^ 4 ) = — u and e^ = e^,i = 1,2,3. Then we can project (fl~6l) onto the observer's 
frame: 

~(a) = (a) . _ (a) _ _ n(a) (i ? A 

z - e dr 2 M dr 2 ~ ( 4 )(«( 4 ) 1 ' J 



1 We will use the indices {1,2,3,4} for basis elements instead of {0,1,2,3} adopted in [3], i.e., 
changing to 4, and changing the ordering. 
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with = Z = e^Z M and -R( a )(4)(/3)(4) = e^u u ej^u 5 R^s ■ Now we introduce the 
second null basis {e„} = {e|, e 2 , e 3 , e^} = {m, rh, 1, k} associated with the observer's 
frame: 

f J ^ (18) 

1 = (u - e (3) ) , k = — (u + e (3) ) . 

One can check that all derivations in Section II of [3] also hold for our non-holonomic 
basis {e a }. Here we only quote those results relevant to our purpose. To begin with, 
all test particles should be synchronized so that Z^ = (they always stay in the same 
local hyper surf ace). For type N spacetimes, the rest of (TTTj) reads 

z^ = (f -a) zw + a x z^\ 

&Q = (± + A + ^Z® + A x Zl 1 \ (19) 
3 

with wave amplitudes of the two polarization modes given by 

A+ = -KeVt, A x = ^m*4, *4 = C a ^ s l a rh^Frh 5 . (20) 

Here ^4 is calculated in the second null basis (fl8"|) . hence different from, but related to 
^4 in (II 01) . Assuming that the observer's frame {e( Q )} is parallel-transported along u, 
i.e., De^/dr = 0, then we have = D 2 (e^ ■ Z)/dr 2 = d 2 Z^/dr 2 , which makes 
( JT9l) easier to solve. 

5. Parallel-transported frames 

Given a radiative spacetime with a principal null direction k and an observer's four- 
velocity u, we can construct the observer's frame {e^)} according to ffl"8l) together with 
the following proposition . 

Proposition 1. [3] Let u be the observer's four- velocity and k be the null vector 
(principal null directions) that satisfy k • u = — ^ . Then there is a unique spacelike 
vector e(3) = A/2k — u. Another null vector 1 is given by 1 = \^2u — k such that 
1 ■ k = — 1. The only remaining freedom are rotations in the transverse plane (e(i), e( 2 )) 
perpendicular to e( 3 ). 

With the null basis {ei, e 2 , e 3 , e 4 }, the metric ([1]) takes the very simple form 

/ 1 \ 

1000 rm 
K 1 1 
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We take the vector field k to be aligned with the quadruple principal null direction e 4 , 
i.e., 

F = (0,0,0,fc 4 ), 

and recall that w M = (it 1 , u 2 , w 3 , w 4 ), u 1 = u 2 in the basis {e a }. The interpretation null 
basis described in (fT8"]l and Proposition 1 has the form 

«n u / u 2 

0,-1,0,- , m^= -1,0,0,- 

IT/ \ If 3 

F = (V2u\ V2u 2 , V2u 3 , V2u 4 + -i— j , (22) 
P = (0,0,0 



v/2u3 



which is unique up to rotations m — >■ e l9 m and trivial reflections. The corresponding 
orthonormal frame in the basis {e a } is given by 



p 

(1) ~ 


V2\ 


p 

(2) - 






P 

(3) ~ 


-(« 


P 

(4) ~ 





U 1 + M 2 

-1,-1,0, — 

It 

u 1 — u' 
-i,i,0, . 

lif 3 



it 3 



(23) 



One can check that these expressions are consistent with Eqn. (18) and (19) of [3] when 
applied to (j2ip with matrix elements rearranged accordingly. 

In general, the frames (e( Q )} and {e^} cannot be parallel-transported along the 
geodesic with u = ea). We can use the following proposition to single out those 
geometrically privileged geodesies along which the interpretation frames are indeed 
parallel-transported. The detailed proof can be found in [8] with no difficulty to be 
generalized for the non-holonomic basis fl6]) with the metric (121j) . 

Proposition 2. [3], [8] Given a geodesic with the tangent vector n M = (it 1 , u 2 , u 3 , u 4 ) 
in the spacetime ( |2T|) . the interpretation null basis ( 122]) and the orthonormal frame ( l23l 
are parallel-transported along the geodesic if 

o = rV (24) 

and 

*ll(r)=irV (25) 
where i?m is the rotation angle for m — > my = e 11? iim in order that Dm^/dr = 0. 
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6. Privileged geodesies 

Now we apply the above results to our twisting type N spacetimes with A. First recall 
that T 2 42 = T 2 43 = T 2 44 = 0. Then the condition (1241) and its complex conjugate are 
tantamount to 

u 1 = u 2 = 0, (26) 

which, by ( TT3|) . leads to 

x = J = 0. (27) 

Along such radial geodesies with fixed x and J, the geodesic equations are quite 
simplified with the first two of f fT4|) given by 

i A R 

= 2AP(d jL y (e_ir + 1)(p/ + 2AJ - 2iCl) ^ (28) 



and its complex conjugate. The rest of (114]) are given by 

° = ii+ (6^i tan 0" 2 ' (29) 

r / 2 A A 2 cosr 2 \ 

° = f + tan 2 {' -3^™-!W7 (30) 

The conserved quantities ( 1T5|) . i.e., first integrals of ( 1291 130|) . read 

= r + u + 4C 2 (0jL) cos 2 -, (31) 

(W + W) or 
= 1 Q r J ii + 4(C , 1 C 2 n - 2C 2 L + C 3 ) cos 2 -, (32) 



in addition to the invariant length ([12]) 

e / A cosr ^ \ . , 

Note that djL and L above are constant for fixed x and J and that is given by (JH]). 

Now we proceed to solve the system fl28H30l) . First note that the equation (1251) 
gives rise to the following two possibilities. 

Case 1. We have u = 0. Then (133|) requires e = 0. This is the case corresponding to 
null geodesies along principle null directions. The geodesic equation (J30l) immediately 
gives us 

2 r 
= f + r tan -, 

which has the general solution 

r = 2 arctan(y4r + B) 
with A, B arbitrary real constants of integration. 

Case 2. Assuming « ^ and a given solution Po(J) to (0J, we have, for ( 1251) to hold, 
= + 2AJ, C x = 0. (34) 
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The first equation above, with its right-hand side being a function of J, shall fix the 
value of J which we call Jo- In fact, combining ( l3~4"j) with we know that Jo is also 
the point at which the second derivative P${J) reaches its maximum value — 10A/3. 
Hence for a given Po(J), we only have limited choices of J for the privileged geodesies 
described in Proposition 2. Nonetheless, the coordinate x can take on any arbitrary 
constant value which we denote clS X — Xq. In what follows, such P = Pq(Jo) and x = x 
will always be assumed. 

To solve the system (129) l30l . one can first solve ( J29|) for r and then combine the 
result with the first integral ( )3~Tj) so as to obtain a third-order ODE for u(t) alone 



U = uu-u 2 -^-ii i -— f u\ A^O. (35) 
3 3oG£ 

Thus for this equation, we obtain in the real domain the following two different waya^] 
to represent its general solution both of which we find suitable for timelike geodesies 
with A > (cf. gl])): 

tanh (l2^ {U ~ El) ) = ±(eAlT+Bl + (36) 

tanh (l^Z {U ~ E2) ) = V 1 + D " tanh ( A2T 2 B2 ) + D2 > (37) 
with Cl = djL(J ) ^ and real integration constants Ai,P>i,Di, Ei (C 2 = AiDi/2Cl) 
and A 2 ,B2,D2,E 2 (C 2 = -A 2 /2C L ^l + D\). By the equation (J2SD, these m(t)'b yield 
respectively 

r(r) = ±2arctan^-^e-( AlT+Bl )(l - Df) - ^ e AlT+Bl ^j , (38) 

Dl 



r(r) = -2 arctan [ " 2 sinh(A 2 r + B 2 ) + D 2 cosh(A 2 r + B 2 ) I . (39) 

One can check that the expressions fl36l 1381) and fl37] [39]) all constitute general solutions 
to ( |29l l30l) in the real domain. In particular when 7J 2 = in ( |37l l39l) . one has a special 
solution 

u(r) = ^(A 2 r + B 2 ) + E 2 , r(r) = 0, (40) 
which is not included in the solution f l36l [381) . From f )33|) . we obtain for both sets of 



3Aj 



~ A ■ ^ 
with e = — 1 for timelike geodesies. In addition, the equation ( 1321) yields trivially 
C 3 = 2CMJ ). 



2 In the complex domain, the expressions (|36|) and (|37|) are in fact equivalent. To see this, one can 
separate a constant from E\ and then use the addition theorem of tanh on the left-hand side of (f3G]) 



and compare it with (|37j) . However in the real domain, this transformation from (|36|) to ([37)1 . for their 
full ranges of real integration constants, cannot be achieved without generally going into the complex 
domain (e.g., note that arctanh(a;) is complex for x > 1). For certain limited ranges of real integration 
constants, the solutions ([36)1 and (|37|) may represent the same solutions. 
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7. Radial timelike geodesies and Killing horizon 



We continue to study timelike geodesies given by ( 1361 1381) and (1371 1391) . For simplicity, 
we only consider the more physically relevant situation with A > 0. Hence from (I4ip . 
one has 

A 1>2 = ±y|, A > (42) 

for timelike geodesies. Note that the range of the hyperbolic function tanh on the left- 
hand sides of (|36|) and ( 137|) is limited to the interval [—1, 1], while the right-hand sides 
are not. Therefore the solution u(r) with its r restricted by the reality condition may 
generally reach infinity at some finite proper time. Nonetheless, whenever this happens 
at the critical value r = r c such that 

e AlTc+Bl +Di = ±1, D\ ^ 1, 

or 

— o ( A 2 T C + B 2 \ 

D|tanhf^-^ -)+D 2 = ±l, D 2 ^Q, 

it always corresponds to 

r(r c ) = 2 arctan( T l) = =f|, (43) 

which are in fact very special hypersurfaces in the spacetime. We can see their 
significance from the metric (pQ) with constant x and J 

■ du fdr H — — cos r du J , — tt < r < tt, (44) 



2C L cos 2 § V 6C 



which indicates that the Killing vector d u is timelike for —tc/2 < r < tt /2 but spacelike 
for — tt < r < —tt/2 and tt/2 < r < tt, and null at r = ±7r/2. Thus that ii(t) diverges at 
r = ±7r/2 under the metric ( 1441) is very similar to that of the Schwarzschild solution in 
the Eddington-Finkelstein coordinates near r = 2m, the fact of which suggests a Killing 
horizoi^ at r = ±7r/2. Inside the region r = ±7r/2, the spacetime is stationary with u 
being a time coordinate (e.g., ( 1401) ). Particularly in the weak field limit \l/4 — > with 
finite A > pQ, one can expect r = ±7r/2 approaching the cosmological horizon of the 
de Sitter universe. 



8. Wave amplitudes 

Now we calculate the wave amplitudes A+ and A x in ( JT9l) . The null bases {e Q } and 
{e,j} (cf. 06]) and (|T8|) ) are related by the Lorentz transformation 

e 4 = A e^, 

d = e ie e 2 + fie 4 , (45) 



-e 3 = A" 1 (e § + Be ie e 2 + Be^ei + Bfie 



3 Duc to cross terms like d£du in the metric with the current coordinate system, it is not at all clear 
how to show r = ±7r/2 an actual horizon for non-radial causal geodesies. Hence we will not dwell on 
this issue. 
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with A = —\^2u 3 , B = —\plu\ and 6 = n. The Weyl scalar \l/ 4 transforms as [5] 

^ 4 = = ( U 3 )2 e ir/2 cos 3 r ^ (4g) 

9 2 

where u 3 = —R(u + 2Lx)/(PdjL) from (Tl3|) (also, use the re-scaled function g(w) from 
(J5]) to see the proportionality to A). Hence we know 

V T 

A x = — u 2 sin - cos -. 
18C£ 2 2 

Substituting either (|36l 138}) or (I3T 1 [39]) for w(t) and r(r) above with (|42|) . we find that 
as long as the limit u(+oo) exists (e.g., when A\ < in (1361) ) and r(+oo) = ±7r, the 
amplitudes behave like 

A + ~ Aexp I -4;/^ r I , ^4. x ~ Aexp I -3W^f r I . (48) 





as the proper time r — > +oo. This means that the gravitational waves are decaying 
exponentially fast with the spacetime locally approaching the de-Sitter universe. 
Hence the cosmic no-hair conjecture [10] is demonstrated with these radial geodesies. 
Furthermore, due to the proportionality to A, the polarization modes A + and A x cannot 
be generally separated from the isotropic background represented by A/3 in ( I19j) . when 
it comes to consider their local effects. 



9. Concluding remarks 

Like their non-twisting counterparts [3], the twisting spacetimes described by (tUHJ) bear 
a similar local interpretation as exact transverse gravitational waves in the de-Sitter 
universe. However, the essential nonzero requirement of the cosmological constant A 
in our radiative solutions (wave amplitudes proportional to A; for A = 0, spacetimes 
becoming flat), to our best knowledge, has not been seen in any other radiative exact 
solution (e.g., Kundt class and Robinson- Trautman class). This distinctive feature 
of these spacetimes suggests that the gravitational waves they represent may not be 
considered as being generated by usual astronomical sources which always radiate 
regardless of the background A being zero or not. Moreover, from the linear theory of 
gravitational waves with cosmological constant, one can indeed identify extra freedom in 
the forms of approximate solutions caused by the A parameter [IT], which is consistent 
with our observation from the exact theory. In fact, according to [UJ, this extra freedom 
in the wave solutions has been described as being "coupling to matter sources with a 
strength proportional to the cosmological constant itself", though the actual physics 
behind this ad hoc coupling is quite unknown. Altogether, both theories suggest a 
quite different role that the cosmological constant (or perhaps, dark energy) may play 
in the process of gravitational radiations, other than simply the "inactive" de-Sitter 
background as it is generally thought of. Since in the inflationary epoch, the effect of 
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A on the cosmology was much greater than at present, we expect that these radiative 
solutions might be of more relevance to primordial gravitational waves from the Big 
Bang. 
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